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We study superconducting states of doped inversion- symmetric Weyl semimetals. Specifically, 
we consider a lattice model realizing a Weyl semimetal with an inversion symmetry and study the 
superconducting instability in the presence of a short-ranged attractive interaction. With a phonon- 
mediated attractive interaction, we find two competing states: a fully gapped finite-momentum 
(FFLO) pairing state and a nodal even-parity pairing state. We show that, in a BCS-type approxi- 
mation, the finite-momentum pairing state is energetically favored over the usual even-parity paired 
state and is robust against weak disorder. Though energetically unfavorable, the even-parity pairing 
state provides an electronic analogue of the 3 He- A phase in that the nodes of the even-parity state 
carry non-trivial winding numbers and therefore support a surface flat band. We briefly discuss 
other possible superconducting states that may be realized in Weyl semimetals. 



I. INTRODUCTION 

Recent progress in our understanding of non- 
interacting Bloch electrons^^ reveals a large class of 
gapped topological phases, t he so -called topological in- 
sulators and superconductor^^. For example, a time- 
reversal-symmetric topological insulator is a band insu- 
lator that cannot be continuously tuned into a trivial 
atomic insulator, as long as time-reversal symmetry is 
respected. A topological insulator is featured by a single 
Dirac cone in its surface state spectrum. Typically, these 
topological phases are realized in systems with strong 
spin-orbit coupling. It is known that when topologi- 
cal insulators are combined with superconductivity via 
the proximity effect^ or via phonon-mediated attractive 
interaction^, the interesting interplay between electron 
pairing and the spin-orbit coupling results in exotic su- 
perconductivity. For instance, when a topological insula- 
tor is doped and turned into a superconductor, an odd- 
parity topological superconductor is obtained^ with a 
single gapless Majorana surface state, which is protected 
by time-reversal symmetry. 

Another type of gapless 'topological matter', the Weyl 
semimetal, is currently being studied intensively and is 
proposed to be realized in experiments^^. Its elec- 
tronic structure has an even number of Weyl nodes - two 
cylyndrical 3D cones that touch at their apex, the Weyl 
point - which carry non-trivial winding numbers ensur- 
ing their stability. These Weyl nodes can be thought 
of as 3D analogs of the two component Dirac fermions 
in graphene and at the surface of a 3D topological in- 
sulators. They exhibit spin- momentum locking and thus 
require strong spin-orbit coupling to be realized. In anal- 
ogy with superconductivity in topological insulators, it 
is natural to expect interesting superconducting states 
to emerge in these systems upon doping, resulting from 
their non-trivial topological winding numbers. To real- 
ize a Weyl semimetal phase requires either time-reversal 
symmetr}fl3] or inversion symmetr}^ to be broken. In 
this paper we concentrate on the inversion-symmetric 



case, in which the two nodes connected by the inver- 
sion symmetry carry opposite chirality. Upon slight dop- 
ing there are at least two disconnected components to 
the Fermi surface around the nodes, shifted in momen- 
tum space from the inversion- symmetric high-symmetry 
points such as the T-point. We will show that the in- 
terplay between the finite-momentum displacement and 
the non-trivial winding numbers around each Weyl node 
leads to interesting superconducting states. 

The finite momentum shift of the Fermi surface mo- 
tivates the study of finite-momentum pairing states 
or Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states^EU. 
FFLO states break translational symmetry and have in- 
teresting physical properties^HUl i n the Weyl semimet- 
als, the center of momentum of the FFLO pairs is fixed 
by the momentum of the Weyl nodes. Similarly, the 
non-trivial winding around the nodes and the broken 
time-reversal symmetry suggests the possibility of realiz- 
ing even/odd-parity BCS states that are electronic ana- 
logues of the 3 He-A phas d 2Q | 21 4 Since the 3 He-A phase 
has nodes with non-trivial winding number which guar- 
antees the existence of a dispersionless surface states 21 , 
the Weyl semimetal in these phases is also expected to 
support zero-energy surface flat bands, similar to a Weyl 
semimetal in proximity to a superconductor 2 -^ 2 ^. 

Surprisingly, when the attractive interaction is com- 
pletely local in real space and represents a phonon- 
mediated interaction, we find from a self-consistent 
mean- field calculation that the fully- gapped finite- 
momentum pairing is energetically favored over the even- 
parity BCS state (both pairing states can be thought of 
as spin- 'singlet' pairings though 'singlet' is not a very 
exact terminology since spin-rotational symmetry is bro- 
ken) and is stable against weak disorder. Hence, there 
is a good chance of experimentally observing these ex- 
otic phases. To be concrete, we concentrate on a spe- 
cific lattice model realizing an inversion-symmetric Weyl 
semimetal and solve the gap equation of the lattice model 
in the BCS approximation. We also discuss the ap- 
plicability of our result to other models realizing Weyl 
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semimetals. 

The proximity effect of an s-wave superconductor on an 
undoped Weyl semimetal has been studied by Meng and 
Balents in Ref. |22j In contrast to this work, where the 
superconductivity is extrinsic, we are interested in the in- 
trinsic superconductivity of the doped Weyl semimetal. 
Note that if the Weyl semimetal is undoped, the intrinsic 
superconducting gap and the critical temperature are ex- 
pected to be vanishingly small since the density of states 
goes to zero at the Weyl point. 

II. MODEL 

The model we consider in this work is given by the 
Hamiltonian 

H = H + V ee (1) 

where V ee is an electron-electron interaction term to be 
specified below. For the kinetic term, Hq, we take the 
minimal two-band lattice modeP^ 

Hq = t(cr x sin k x + cr y sin k y ) + t z (cos k z — cos Q)cr z 

+ m(2 — cos k x — cos k y )a z — fi. (2) 

This model realizes a Weyl semimetal with two Weyl 
points at momenta P± = (0,0, d=Q). a x,y,z are the Pauli 
sigma matrices (for later use we define <r° to be the 2x2 
unit matrix), t and t z siriQ are the Fermi velocities at 
the Weyl points in the x, y and z directions respectively. 
Without a loss of generality, we assume t = t z sin(Q) 
such that the Fermi velocity around the Weyl points is 
isotropic. We have explicitly included the chemical po- 
tential fi in the kinetic term. We are primarily inter- 
ested in the parameter range < |/z/t| <C Q, when the 
Fermi surface consists of two disconnected spherical com- 
ponents around the Weyl points (see Fig. [I]). In this case, 
the states on the Fermi surface have spin-momentum 
locking similar to the surface states of a strong topo- 
logical insulator. This property will play an important 
role in our discussion of pairing states below. 

The electron-electron interaction is short ranged and 
takes the form 

V ee = Vq ^ riiTli + V\ ^ n i n J = V ^ n k U -k^ ( 3 ) 

i (ij) k 

where rii = J2 a c l a c i,o- 1S the number of electrons on 
site z, and the second sum is over nearest neighbors only. 
V(k) = V + Vi(cosfc x + cos ky + cosk z ) is the Fourier 
transform of the real- space interaction. Vb represents 
a phonon-mediated attractive on-site interaction and V\ 
the nearest-neighbor interaction. We are mainly inter- 
ested in the case Vo < and \Vo\ ^> |Vi| when electrons 
form Cooper pairs and condense. The phenomenologi- 
cal interaction term (|3| captures the tendency of d-wave 
pairing for V\ < and Vo > in the context of high-Tc 
superconductors such as cuprates^. 



The point group symmetry of the model Hamiltonian, 
H, is C 4h = {l^Cf 1 77/ = 0, 1; 774 = 0, 1, 2, 3}, where 

I :a z H(-k)a z = H(k), 
C 4 :S^H[R n/2 (k)]S = H(k) 1 (4) 

with S = ^( cr ° + ^ (jZ ) an d R-k/2 a rotation by an angle 
7r/2 around the z-axis. / is the inversion symmetry that 
takes r — >• —r. Each spatial rotation is accompanied by 
an equal spin rotation, manifesting the spin-momentum 
locking due to the presence of strong spin-orbit coupling. 
The CVrotation symmetry in the xy-plane is not required 
to realize a Weyl semimetal, but is present in this model 
and similar lattice rotations are present in other models 
we discuss later. 

To make a connection to previous worlP^^ and to ob- 
tain a general understanding of Weyl semimetal phases, 
we derive the low energy effective theory corresponding 
to the lattice model jl]). Expanding H(k) in the small 
momentum q = k — P± around the two Weyl Points de- 
noted by ±, we obtain 

h = £ c t(fc)i/ (fc) c (fc) » ^ a (q)ha(q)Mq). (5) 

k a=± 

The effective kinetic Hamiltonian h±(q) is given by 

h±(q) = t(q x a x + q y a y =F q z a z ) - /i. (6) 
Similarly, we obtain for the interaction term 

V ee = V ab ' cd (q)^l a (k + «)< T (P - q)^ c ,r(P)Mk), 
k,p,q 

(7) 

where roman letters denote the nodal indices ± and <j, r 
are spin indices. Here and henceforth, repeated indices 
are summed over. In the BCS channel (see App. [A] for 
details) 

v ee = g y a ^ c ^t )<T( fc )V) t >r( _fc^ C)T( _ r)V;(i)(7(r)) (8) 

kX 

with 

v -+;+- = v+-,-+ = Vo + 3Vl _ _ f)2 ; 

= Vo + 2V, + V ± + Vj+ (9) 
V+--+- =V + 2V X + V± + Vf, 

and 

V ± = -^(k ± -l ± ) 2 , (10) 

V+/V! = [1 - X -(k z - l z f\ cos 2Q + (k z - l z ) sin2Q, 

Vf/V! = [1 - \{k z - l z ?}cos2Q - (k z - l z )sm2Q, 

where k± = (k x , k y ,0). These expressions will be used in 
the next section. 
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FIG. 1. A schematic diagram of the spin texture around the 
Weyl nodes in momentum space and the pairing states, (a) 
The spin direction of the eigenstates is given by thick arrows. 
The double-headed arrows labeled with (1) and (2) indicate 
the partner states in the BCS pairing. The spin state is maxi- 
mally anti-parallel for (1) and parallel for (2), indicating that 
there will be nodes in the latter case if the pairing is in the sin- 
glet channel. Contrary to the BCS pairing, the FFLO pairing 
(3) connects two states within the same node ('intra- nodal' 
pairing) . The two states connected by the FFLO pairing have 
the opposite spin directions, (b) Position of the nodes for the 
even-parity state. The nodes of the same chirality are on 
the same component of the Fermi surface, with their partner 
nodes of opposite chirality on the other. The filled circle rep- 
resents a node of chirality +1 and the crossed circle represents 
a node of chirality —1. Hence, there are four nodal points on 
the Fermi surface of the even-parity paired BCS state. 



III. MEAN FIELD THEORY AND PAIRING 
CHANNELS 

We treat the interaction term V ee in a mean-field ap- 
proximation and solve the resulting gap equations self- 
consistently. In addition to the more standard BCS par- 
ing, we also study finite-momentum or FFLO pairing. 





IRR 




pairing function 


r 1 


A 


1 


Z 2 (ia y ), (X - iY)(a° + a z ), (X + iY)(a° - a z ) 




B 


-1 


XY(ia y ), (X + iY)(a° + a z ), (X - iY)(a° - a z ) 




E + 


i 


(X + iY)Z(ia y ), Z(a° + a z ),XYZ(a° - a z ) 


r 3 - 


E- 


—i 


(X - iY)Z(ia y ),XYZ(a° + a z ), Z(a° - a z ) 



TABLE I. Symmetry classification of the BCS pairing order 
parameters for the model |2]) according to the different irre- 
ducible representations (IRRs) of the group C^h- X, Y, and 
Z are basis functions for the momentum-space pairing func- 
tion which we take to be sinp^, sinp y , and sinp z respectively, 
and can be realized by nearest-neighbor pairings. Among the 
pairing functions, ia y means singlet pairing, a x is the spin- 
full triplet pairing, cr° + a z the triplet pairing for polarized t 
spins and cr° —cr z the triplet pairing for polarized \. spins. The 
paired state Y 2 has nodal lines, while the other three states 
are gap less. 



The standard BCS pairing term connects two Weyl nodes 
in the effective theory. The explicit f orm of and A ar 
can be found in Table [I] and in Eq. (A6) in App. |A} The 
self-consistent gap equation takes the form 



Af T (p) = J2v ab ' cd (p- fe)(V c ,r(-A)^(fe)), (12) 



where the expectation value is taken with respect to the 
mean-field superconducting state (see App. [A] for more 
explicit expressions for the gap equations). 



B. FFLO pairing 

In the doped Weyl semimetal, the Fermi surface is 
formed around the Weyl points P±, and it is natural to 
expect a finite-momentum pairing to compete with the 
standard BCS-paired states. We therefore introduce a 
FFLO state with a center of momentum at 2P±, which 
paring function satisfies 



Ap FLO (r) cx exp(2zP + • r) ± exp(2zP_ • r) 



(13) 



A. BCS pairing 

The symmetry classification of different BCS pairing 
order parameters in a doped Weyl semimetal, according 
to the lattice symmetry Q, is summarized in Table [I] (see 
App. [C] for more detailsjT There are three fully-gapped 
BCS pairing order parameters (T 1 and T 3j± ) and one that 
has nodal lines (T 2 ). 

In the continuum theory, the pairing terms of Table [I] 
take the form 



y> ffT (£)4(fc) c t(-fc) 



£A^«1VC(9)< T (-9), 
q 

(ii) 



The self-consistent equations for these pairing order pa- 
rameters take the form 

A ar (p;±P) = J2V(P- fe)(^±,r(-*)^±,a(fe)), (14) 



where the two nodes ± are decoupled. These FFLO 
states correspond to the intra-node pairing, in contrast 
to the BCS case which is inter-node pairing (see Fig. T]). 
The two states of the pairings Ap FLO in Eq. ([l3| with a 
relative phase of ±1 between the two components of the 
Fermi surface have the same mean-field energy since the 
two nodes are decoupled in the mean-field theory. 
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c 4 


Pairing function 


r 1 


1 


Z 2 (ia y ), (X - iY)(a° + a z ), (X + iY)(a° - a z ) 




-l 


XY(ia y ), (X + iY)(a° + a z ), (X - iY)(a° - a z ) 


p3,+ 


i 


(X + iY)Z(ia y ), Z(a° + a z ),XYZ(a° - a z ) 


r 3 - 


—i 


(X - iY)Z(ia y ),XYZ(a° + a z ), Z(a° - a z ) 



with a gap equation 



TABLE II. Classification of the FFLO states of superconduct- 
ing Weyl fermions based on the lattice symmetry C4. The 
notation is the same as in Table [J We assume that the cen- 
ter of momentum for the pairing is at P±. Note that the 
symmetry is only C4 on the xy-pl&ne without the inversion 
because 'inversion' is already encoded by the ansatz Eq. (13). 



This classification is essentially the same as that of BCS-type 
pairing order parameters. 



C. Mean field energy 

Having identified the possible superconducting states, 
we compute their free energy by solving the self- 
consistent gap equations numerically. We are interested 
in the case \Vq\ ^> |Vi| with Vo < where the spin sin- 
glet is preferred. We denote the pairing term oc i<j y 
in Table [i] as 'singlet' and the other terms ex iaa y as 
'triplet'. The singlet and triplet components have a dif- 
ferent dependence on the interaction parameters Vo and 
V\. The gap of the triplet components depends only on 
the value of V\ , while the singlet component depends only 
onFo+3Fi « 1 K ) orF + ^i(5+cos2Q)/2 « V . We there- 
fore consider in the following only the singlet component 
oc icjy of T 1 for the BCS and FFLO states. 

For these two states the BCS mean-field approximation 

is 

H = H + K P e air , (15) 

where Ho is given by Eq. ([2| and V^ air is the effective 
projected pair potential derived from the lattice interac- 
tion in Eq. ([3]). For the T 1 -BCS state we have 



k,p 

pj.=^(k) T x ia v ip*(-k), 



(16) 



and the gap equation is 
A Uses 



> a , a (£0(rTV^T / W-fc)>- (17) 



4 J k 

For the r 1 -FFLO state we obtain 

ypair = _ f/FFLO Y-ptp. 

k,p 

C/fflo = Vo + Wi, 



(18) 



FLO 



(^aAk)(-^ y ) a ^aA-k)) (19) 



In this standard BCS-type approximation, we can eval- 
uate the energy E of the pairing states with the pairing 
amplitude A^ = — U (P_%) (with the effective pairing in- 
teraction strength U) by computing 

E = E e \ + E sc , 
E e \= ^2 e sc (k)n e [e sc (k)} 

kedQ,e sc <0 

^2 e fe (fc)n e [€ fe (fc)], 



E*. 



kedQ,e is <0 

: - £ -4^ + h.c. 



2U 



(20) 



kedn 



Here n e is the filling of the electron for the state at k 
of energy e(k), e sc is the energy of the filled band of the 
BdG quasiparticle with mean-field gap A, and ef s is the 
energy of the filled bands of the free Weyl electrons with- 
out pairing, i.e. the energy of the normal state. Thus, 
the second line of Eq. ( 20 ) represents the energy gain of 



the superconducting state relative to the normal state by 
opening up a gap near the Fermi surface. The last line of 
Eq. ( 20 ) represents the contribution from the pairing in- 



teraction labeled by the momentum k. The range of the 
summation is restricted to a shell dtt around the Fermi 
surface, which width is determined by the strength of the 
attractive interaction. 

In Fig. [2] we plot the mean-field energy for the two 
parings, as obtained from Eq. (20), as a function of 
the interaction strength Vo. The L^-FFLO state has a 
larger gap than the T 1 -BCS state and is energetically 
favored. This result can be understood by consider- 
ing the spin-momentum locking around the Fermi sur- 
face. For the even-parity pairing state, the state |fc,a) 
(a is the spin state) is paired with the inversion partner 
state I — k,<j z a). The pairing amplitude is of the form 
~ (c^ (k)ia y c* (— k)) which takes the maximum value if 
the two states at k and —k have opposite spins. How- 
ever, the spins at k and —k are not anti-parallel and even 
become parallel at the poles (which is the origin of the 
nodes, see Fig. [I]) which tends to reduce the supercon- 
ducting gap. In contrast, the FFLO state connects the 
states \k + Q,a) and | — k + Q,/3) via the spin singlet 
channel with f3 = —a (anti-parallel spins). A gap opens 
up everywhere at the Fermi surface with a larger gap 
than the even-parity BCS state. This is very similar to 
the surface of topological insulator which we discuss in 
Appendix B. 
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FIG. 2. Mean-field energy E of the even-parity I^-BCS and 
the r 1 - FFLO states as a function of interaction strength 
Vo. Other model parameters used to obtain this plot are 
fji/t = 0.3, Q = 0.7, Vi = 0, and dtt = 0.2 



IV. DISCUSSION 

In this section, we discuss the nodal structure of the 
r 1 -BCS state and the effect of the disorder on the T 1 - 
FFLO state. Like in the last section, we consider only the 
singlet components of these states, assuming |Vq| ^> \Vi\. 



A. I^-BCS state 

The pair potential term in the mean field Hamiltonian 
of the singlet component of the T 1 -BCS state is 

H*£? = J2 Act (fc)(i^)^4(-fc) + h.c. (21) 

k 

= j2^lM(T x ) ab (i° v r^U(- ( fi + h - c - 



The second form is obtained in the low energy theory. 
This superconducting state is an even-parity state and 
has four point nodes on the Fermi surface at q x = q y = 



and q z = A 2 + ji 2 (see Fig. The gap remains 
closed even when the triplet pairings of the I^-BCS state 
in Table [J are included. 

The two nodal points near Weyl node P + (P-) carry 
a winding number of +1 (-1). To demonstrate this we 
write down the Bogoliubov-de Gennes (BdG) Hamilto- 
nian H = Y^k^k^k^k for ®k = ( c k' icTl /c lfc) 7 • In the 
continuum limit at P + (similar expressions are obtained 
for P_) 



H 



h + (q) Aa° 
.Act -h_{q) j 



(22) 



with h±(q) defined in Eq. The quasiparticle spec- 

trum corresponding to this BdG Hamiltonian is 

E(q) = ±[q 2 + A 2 + /i 2 ± 2(A 2 q 2 z + ^q 2 ) 1 ^ 2 (23) 



which has nodes at q x = q y = 0, q z = ±^ A 2 + /i 2 , 
with chirality of +1. Near the nodes Ig^l, \q y \ <C \q z 
we obtain the anisotropic Weyl spectrum 



E(q) « ± 



(g z ±V^T7?) 2 + 9i(l 



A 2 



both 

1/2 

(24) 

with q± = (q x ,q y )' At zero chemical potential, this is 
similar to the results of Meng and Balents 22 who con- 
sidered the proximity effect of undoped Weyl semimet- 
als. The effect of nonzero chemical potential is to simply 
shift the Weyl nodes located at q z = ±A at \i = 0, to 
Qz 



Because of the non-trivial winding number carried by 
the nodes, the nodal points are robust against small per- 
turbations. The only way to gap out the nodes is to 
undergo a pair-annihilation of nodes with the opposite 
winding numbers, and the nodal points are topologically 
stable as long as they are separated enough in momentum 
space. Strikingly, this nodal structure implies that there 
will be a zero-energy state on the surface which should be 
detectable in experiment. This is similar to 3 He- A which 
is an odd-parity pairing state, while our superconducting 
phase is realized by the even-parity pairing. 



B. T^FFLO state 

The singlet component of the T 1 - FFLO state is fully 
gapped with a mean-field pair potential 



if™ = Act (k + P + )(iay)^cl(-k + P + )±(P + o P_) 

(25) 

with center-of- momentum of 2P±. In the low-energy 
theory, it can be represented by the intra-node pairing 

It is known that some two-dimensional FFLO states 
with strong spin-orbit coupling and paral lel magnetic 
field are unstable against weak disordei^ 2 ^. In contrast, 
the FFLO state discussed in this paper is found to be ro- 
bust against weak disorder. In fact, the structure of the 
FFLO state Eq. (13) and Eq. (25) is more similar to the 



even/odd-parity state of the doped topological insulators 
studied in Ref. 27 than usual FFLO states in the two spa- 
tial dimension. This similarity is manifested if we write 
down the pairing for the Weyl fermions in the continuum 
limit in the helicity eigenstates 

A± oc e^ + m + {-q)) ± ftM$)V- ("$)>], (26) 

which corresponds to Eq. (5) of Ref. 27, Within this 
Cooper channel, we add a scalar disorder term to the 
Hamiltonian 



imp 



Vi 



imp ^ 

k,peFS 



"1 



E 

q,leFS 



^Vl,Af (27) 
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The matrix element V"i is given by 

q,l ° J 



yab . 

g,i 



V ir 



W) (0) 

Mb <#>, 



(28) 



where we have used the standard normalized spin state 
q • a\q) = \q) and q • S\q) = \q) with q = q/\q\ and 
q = — q z )/\q\- With this impurity scattering, the 
self-energy can be worked out in the self-consistent Born 
approximation, and we find that the correction to the 
self-energy and the Cooperon diagram due to disorder 
are exactly of the same form as obtained by Michaeli and 
FiP3. In fact, the only difference between the FFLO state 
A± in Eq. (25) and the even/odd-parity paired states of 



Ref. [27] is phase factors in the matrix elements of V^j 
in Eq. (28) which does not show up in the corrections 



to the self-energy, the Cooperon diagram, and the pair- 
ing susceptibilities. Thus we conclude that the critical 
temperature of T 1 - FFLO states is not affected by the 
disorder and thus T 1 - FFLO state is robust. 



V. CONCLUSION 

In conclusion, we have studied the possible super- 
conducting states of doped inversion-symmetric Weyl 
semimetals. We considered a concrete lattice model real- 
izing a Weyl semimetal and found that the FFLO state 
has a lower energy than the even-parity state if the in- 
teraction is phonon-mediated, and the phase is argued to 
be stable against disorder. Though the even-parity state 
is less favored in energy than the FFLO state, it inter- 
estingly provides an electronic analogue of 3 He- A phase. 

We remark briefly on the implication of our work for 
superconducting states of Weyl semimetal models other 
than the one studied in this paper. Among the many pro- 
posals for the Weyl semimetal phase, we restrict ours elves 
to the models based on the topological insulator d 12 * 14 1 
with the time-reversal breaking perturbation. 



H vr z a • k± + r x k z 



ma 



The typical symmetry of the model is Ix C n ( x M, Mirror 
symmetry) where / is the inversion symmetry and C n is 
the n-fold lattice rotation symmetry along a certain axis 
(for the model based on Bi2Se3 14 , we have n = 3). Due 
to the strong spin-orbit interaction, the spatial symme- 
try operation involves the spin/orbital operations, e.g., 
I : k — >► —k should involve r y , and r y H{-k)r y = H(k) 
(the lattice rotation will involve a spin rotation). Note 
that these symmetry considerations already manifest the 
similarity between the realistic model and the simplified 
model Eq. and this similarity becomes much clearer 
if we go to the low-energy theory of Eq. (29). It is not 



difficult to confirm that the low-energy theory is identical 
to Eq. (J6|, and hence we will have similar superconduct- 
ing states, FFLO and electronic analogues of 3 He-A, in 



the more realistic model. Hence, we predict that the su- 
perconducting states we found should show up in other 
proposals for Weyl semimetals. 



Note that FFLO state Eq. (25) shows a density modu- 



lation pinned by the momentum of the Weyl nodes (which 
is reminiscent of the field-induced charge density wave^ 
of the Weyl semimetals). In the presence of multiple Weyl 
nodes (for example, the iridates have 24 Weyl nodes^ 
in momentum space), there will be multiple centers of 
momenta for the FFLO state and this directly implies 
that there will be interesting density modulation patterns 
which are fully determined by the position of the Weyl 
nodes. (This is true at least at the level of mean-field 
theory which ignores the effect of 0(A 4 ) terms in the 
Landau-Ginzburg theory. 0(A 4 ) terms can potentially 
melt this pattern). 
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Appendix A: Continuum Weyl fermions: 
interaction, and gap equation 

In this appendix, we derive the equations related to 
the continuum theory from the microscopic lattice model 
Eq. ([5]). First of all, let us derive the interaction Eq. ([9| 
for the Weyl fermions. The interaction (|3| in the mo- 
mentum space is 

H = E V k c t(k+P)4(q- k)c T (q)c a (p), (Al) 

k,p,q 

(29) with Vg = V + V 1 (cos(k x ) + cos(k y ) + cos(fc^)). For the 
BCS pairings where the center of momentum is at zero 
q — — p, we obtain the pairing potential 



H = 



E 

p,q 



V p ^t(p)cH-p)c T (-q)c a (q). (A2) 



As the electron operators that we are concerning are lo- 
calized near the Weyl points, we expand the electron op- 
erators near the Weyl points. This can be easily done by 
plug ging p = k±P and q = l±P with P = (0, 0, Q) into 
Eq. (A2). For example, the interaction (|A2|) includes the 



interaction 



~>V, 



(fc+P)-(r+P)^+,a( fc )^, r (-fe)^-,r(-0^+,a(0 



-(k - Z)^, <r (*)^L, T (-*)^-,r(-Z>+ l a(0 > 



(A3) 



7 



which allows us to identify V + ' + (k — I) = 
^(£ +j P)_(f + P)- Similarly, we can identify V (/c — 

= ^(£_p)_(f_p)> ^ + ' +(k - I) = V(fc_p)_(r +j p) 5 and 
F H ' +- (/c — /) = V^ + ^_^_^. After this identification, 
it is straightforward to expand for the small /c, / to obtain 

v -+.+- = v+ -,-+ = Vo + Wi _^±^_ f )2; 

= y + 2 yi + v ± + V+, 
V+-'+- = V + 2V 1 + V± + Vf, 

v ± = -^{k ± -T ± )*, 

V+ = ^(cos(2Q)(l ~\{k z - l z f) + (k z - l z ) sin(2Q)), 

Vj," = Vi(cos(2Q)(l - l -{k z - l z f) - (k z - l z )sm(2Q)). 

(A4) 

Next, we discuss the mean-field pairing terms in the con- 
tinuum theory derived from the lattice model. The pos- 
sible pairing states from the interaction Eq. ([3| are listed 
in the table [IJ The typical form of the pairing term can 
be represented by 



H pa[l = ^Ac^k)T(k)c^-k) 
= E^cUk)T ar (k)cU-k). 



(A5) 



(that is, A aT (k) = AT ar (k)). As we did for the interac- 
tion terms, we expand the pairing terms near the Weyl 
points by 

ct(k)T aT (k)cU-k) 

= ^>ir+-(+,p)< r (-p) 

T 2 = (p x + iPy)T X (<T° + d Z ), (p x " iPy)T X (a° - (J*)' 

r 3)+ = (p x + ip y )Tyay,iTy(a° + a z ); 

r 3l - = (p* - ip y )r y ay, ir»(a° - a z ). (A6) 

As we now have the explicit form of the pairing term 
and the interaction for the Weyl fermions, it is now 
straightforward to obtain the gap equation for each 
ansatz by solving 

Af T (p) = V ab ' cd (p- kMcA-kWdAk))- (A7) 



iY)(a x — icr y ){icr y ) which can mix each other (see table 
A ' = \ Y.UV X ' ,Cd $) - iV v ' cd {k)][Xcd{k) + iYcd(k)}}, 

k 



where we define the following compact notations 
Y,^ d n(r x )a bPx V^ cd {k-p) 



(A9) 



V x ' cd {k) 



pedQ fx 

b.cdt 



v ^ (l)= ^(r*UP,vy*-r> (A10) 

z^pedn Py 

where dQ is the thin shell around the Fermi surface. The 
width of the shell is determined by the phenomenolog- 
ical parameter (Debye frequency) defining the electron- 
phonon coupling, e.g.,t x d£t is the characteristic energy 
of the phonons. 

X cd (k) = ^cA-k)[-i° v ° x } Pa ipdAk)) 

Y cd (k) = (iP c A-k)i-i(T v v v f a ipdAk)) 

Z cd (k) = (i>cA-k)i-iv v <J z f a i>dAk)) (All) 

with the expectation value taken for the mean-field su- 
perconducting state. 

For T 2 pairing (C4 = — 1), there are two pairing chan- 
nels - A 7 (X + iY)(a x + i(j y )(i(jy) + A H (X - iY)(a x - 
ia y )(ia y ) which can mix each other (see table [i]), and the 
self-consistency requires 

A ' = l T,^''" 1 ^) + iVy' cd (k)}[XAk) + iY cd (k)}}, 

k 

A " = I Y.UV X ' ,Cd $) - iV v ' cd {k)][XAk) ~ iYcd(k)}}, 

(A12) 

For r 3 '^ pairing (J = — 1 and C4 = the analoguous 
expressions can be derived by similar methods. 



Appendix B: FFLO state on Surface of topological 
insulators 



For example, the gap equation fo r th e singlet pairing 
component A^ = AfT 1 )^ in Eq. ([A6| is 



A=lZk(T X )abV ab ^ d (p-k) 

{^ c A-k){-i(J y Y a ^a(k))' 



(A8) 



The similar expressions of the gap equations hold for 
other pairings and the results are following. For the 
triplet pairing component in T 1 (C4 = 1), there are two 
pairing channels ~ Ai(X-iY)(a x +icr y )(icr y ) + A n (X + 



In this appendix, we will propose a possible FFLO 
state from the surface state of the strong topological in- 
sulator under the parallel magnetic field. As far as the 
magnetic field is in-plane, the only coupling from the 
magnetic field to the surface state is Zeeman coupling. 
Then, the Zeeman coupling will shift the Fermi surface 
uniformly along the direction to the magnetic field. The 
low-energy theory is 



H = va ■ (k + gB/v) - \i. 



(Bl) 
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IRR 


c 4 


basis functions of A(k) 


Nodes? 


A u 


1 


Z 


N/A 


B u 


-1 


XYZ or (X 2 - Y 2 )Z 


line nodes 


E u + 


i 


X + iY 


N/A 


E u - 


—i 


X-iY 


N/A 



TABLE III. Symmetry classification of distinct BCS pair- 
ing order parameters in (C3), corresponding to different 



odd-parity irreducible representations (IRRs) of point group 
C^h- Here (X,Y,Z) are basis functions for the momentum- 
space pairing function A(k), denoting e.g. the function 
(sm(p x ),sm(py),sm(p z )) or other functions with the same 
symmetry. 



If there is a phonon-mediated surface transition toward 
superconducting state, the superconducting state should 
be at the finite center-of-momentum pairing. We model 
the phonon-mediated attractive interaction as 



5H = uY j np 



(B2) 



which is local in the real space and uniform in the mo- 
mentum space. This implies that the pairing will be uni- 
form in the momentum space, hence we can single out a 
single pairing state A(r) ex exp(— 2iQ-f) with Q = gB/v, 

H paiT = A(Q)cl(Q + k)(ia^cl(Q -k) + h.c. (B3) 

Hence, we have shown that FFLO state can show up, 
at least in the mean-field theory, in the surface state of 
strong topological insulators under the in-plane magnetic 
field. What would be the effect of weak disorder to this 
phase? By following the discussion in reference [26], we 
conclude that this FFLO state should be robust against 
the weak neutral disorder (this problem corresponds to 
the problem where the disorder scatters electrons only 
within the single Rashba band in the reference). 



Appendix C: Symmetry classification of pairing 
order parameters in two-band model ([5]) 

In this section we classify possible BCS-type pairing or- 
der parameters in two-band model Q with interactions, 
as long as interaction terms do not break the point group 
symmetry of tight-binding model These differ- 
ent pairing order parameters are characterized by distinct 
irreducible representations of group C4 generated by 90 
degree rotation along z-axis. 



When the Weyl semimetal described by two-band 
model ([2| is doped slightly with \i > and \jJ>/t\ « 1, 
the Fermi surface consists two electron pockets around 
Weyl nodes ±P = (0,0,±Q). The tight-binding model 
([2]) can be diagonalized into H(k) = d^-a — fi = 

U l( £ k aZ ~ where U % = 



is a 2 x 2 uni- 



tary matrix. Hence the low-energy degree of freedoms on 
the Fermi surface e% = fJ> are 



(Cl) 



and according to inversion symmetry we have U_r 



U^(j z and 



f-k ~ U k C -k,t V k C -kX 



(C2) 



For a general electronic model with C^h symmetry (but 
no time reversal), the Fermi surface is nondegenerate. 
Therefore a generic BCS-type pairing term is written as 



#pair = £( A (£)/ fc ! /lfc+ h.C.), 



(C3) 



Apparently the pairing order parameter A(— k) = — A(fc) 
is a odd-parity function (instead of a 2 x 2 matrix). One 
can always choose a gauge so that under C4 symmetry 
in (El): 



C 4 



k y , k z ) — > (k yj k Xj k z ) 



the eigenvector (wr, vA transforms as 



Aj C 4 k 



J C 4 k 



(C4) 



(C5) 



In this specific gauge the operator /g transforms trivially 
under rotation C4 and inversion /, hence the order pa- 
rameter A(fc) should form a one-dimensional odd-parity 
(irreducible) representation of the symmetry group C±h- 
All possible different order parameters are classified by 
their symmetry and listed in Table III The B u state 
has nodal lines along which there are gapless excitations, 
while the other three BCS paired states are fully gapped. 

The symmetry classification of BCS pairing order pa- 
rameters are generally true for any system with C^h sym- 
metry but no time reversal (so that the Fermi surface is 
nondegenerate). In the specific two-band model we 
can relate the low-energy eigenband operators /g to origi- 
nal electron operators a through (Cl ), and the distinct 



pairing functions in the original electron basis are sum- 
marized in Table HJ 
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